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Abstract 

The variation of the velocity-derivative skewness of a Navier-Stokes flow 
as the Reynolds number goes toward zero is calculated numerically. The value 
of the skewness, which has been somewhat controversial, is shown to become 


small at low Reynolds numbers. 


We write the skewness S of the velocity derivative Bu^ax! 


as 



where uj is the velocity component in the Xi-direction and the overbars 
indicate averaged values. If the skewness (or the skewness factor) is a 
measure of the nonlinearity of the flow, one might expect it to approach zero 
as the Reynolds number becomes small, since the Navier-Stokes equations should 
be linear at low enough Reynolds numbers. Batchelor and Townsend [l] have, in 
fact, given a relation between the nonlinear term in the vorticity equation 
and the skewness . 

Nevertheless, some plausible analyses have indicated that the skewness 
factor may not approach zero at vanishingly small Reynolds numbers. In 
particular, Refs. [2] and [3], where the shape of the energy spectrum was 
assumed to remain similar with increasing time, concluded that the skewness 
factor approaches a nonzero value in the final period of decay. That value 
did not differ greatly from those at earlier times. On the other hand an 
analysis where the assumption of similarity was not invoked [4 and 5], and an 
experiment [5], indicated that the velocity-derivative skewness approaches 
zero in the final period. 

Here we use a numerical solution for a decaying Navier-Stokes flow to 
determine how the skewness changes with time. The Navier-Stokes equations for 
an incompressible flow can be written as 

3^ a(u t u k ) 1 3^ aV 

3t Sx^ p 3x. + u ax^Sx^’ (2) 
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( 3 ) 


where the pressure is given by the Poisson equation 


1 _ 3 <Vk> 

P 3x 2 9x a 


3x a 3x k ' 


The subscripts can have the values 1, 2, or 3, and a repeated subscript in a 
term indicates a summation, with the subscript successively taking on the 
values 1, 2, and 3. The quantity uj is an instantaneous velocity component, 
x[ is a space coordinate, t is the time, p is the density, i> is the 
kinematic viscosity, and p is the instantaneous pressure. As in [6] the 
initial velocity is given by 

u[ = aj cos qx + bj cos r*x + Cj cos s-x, (4) 

Where ai - k(2,l,l) . bi = k(i ,2 ,1) , ci = k(l , 1,2), 

qi - (-l,l,l)/x 0 , ri = (l,-l,l)/x 0 , sj = (1,1 ,-l)/x 0 , (5) 

k is a quantity that fixes the initial Reynolds number at t = 0, and x 0 is 
one over the magnitude of an initial wavenumber component. The initial 
pressure is not specified since it is calculated from Eq. (3). Equations (4) 
and (5) satisfy continuity, and Eqs . (2) and (3) insure that continuity is 
maintained. The boundary conditions are periodic, with a period of 2ffx 0 . 

It has been shown in [7] that 128 3 grid points appear to resolve our 
flow, even for times when the Reynolds number is relatively high (initial 
Reynolds number - 1,000). Since we are mainly concerned here with the final 
period of decay, where the Reynolds number is low, the results of most interest 
should be very well resolved spatially. As in [6], numerical stability 
limitations force the timewise resolution to also be good. 

Figure 1 shows the calculated time evolution of the velocity-derivative 
skewness. At the earlier times the skewness oscillates, having an average 
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value somewhere around 0.5. That type of behavior was also obtained for the 
numerical calculations in [8]. For later times S goes monotonical ly toward 
zero as the fluctuations decay. 

The calculated variation of the skewness with the microscale Reynolds 

_ 1/2 
2 

number R\ - u X/u is plotted logarithmically in Fig. 2, where 
~2 

u = u jUj/3 and X is the Taylor microscale. Except at the larger Reynolds 
numbers, the skewness goes toward zero monotonical ly as R^ decreases. The 
limiting rate of approach to zero for the numerical solution is about the same 
as that for the experiment in [5] (S R\ 14 ). 

The nonsimilar analyses in [4] and [5] give skewness factors that 
approach zero as the Reynolds number decreases, as do the skewness factors in 
our numerical solution and the experiment in [5]. This is in contrast to 
similarity solutions [2 and 3]. The results suggest that the use of similarity 
assumptions for calculating skewness factors in the approach to the final 
period of decay may not always give realistic trends; the change in shape of 
the energy spectrum with time seems to have an important effect. 

As a final observation we note that at early times the fluctuations in 
our flow are time-dependent and display sensitive dependence on initial 
conditions [7], The flow at early times therefore contains important 
ingredients of turbulence [9]. That is evidently not the case at later times 
in the approach to the final period of decay. In [10] it is shown that, at 
least for forced flow, the fluctuations become time- independent at Reynolds 
numbers ~5 and below. The fluctuations are then completely spatial. 

Time- independent fluctuations also occur for our present flow in the approach 
to the final period of decay, where the Reynolds numbers are on the order of 
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one and less (see, for example, Fig. 3). Thus the flow in the approach to the 
final period cannot be considered truly turbulent. Batchelor and Townsend 
[11] recognized early the nonturbulent character of the flow in the final 
period of decay. Our results indicate that the flow in the approach to the 
final period is also not turbulent in a strict sense, since the fluctuations 
are t i me- independent , being only spatial. The presence of spatial 
fluctuations in the absence of temporal fluctuations is shown by the fact 
that ui is different at the two points considered in Fig. 3, even when 
temporal fluctuations are absent. 

Acknowledgment - I should like to thank Frank Molls for his work on the 
numerical calculations for this study. 
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Figure 1 .—Calculated time-evolution of 
velocity-derivative skewness for a de- 
caying flow. Overbars designate space 
averages. 



Figure 2. — Variation of velocity-derivative skew- 
ness with microscale Reynolds number. 



At grid center, 
x* =x 1 /x 0 =X 2 = X 3 = k 


\ At x * 9rc/8, 


x 2 = 21ic/16, x 3 = 23jc/16 


l. I l l l l l 


8 1.0 1.2 14 1.6 1.8 2.0 

t'= (V/Xq) t 


Figure 3. — Calculated time -evolution of un- 
averaged velocity components for a decay- 
ing flow. Ordinates normalized by initial 
condition. 


8 


IVIASA Report Documentation Page 

National Aeronautics and 

Soace Administration — _ — — 

1 Report No. 2. Government Accession No. 

NASA TM- 103 186 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

Numerical Solution for the Velocity-Derivative Skewness of 
Reynolds-Number Decaying Navier-Stokes Flow 

a Low- 

5. Report Date 
June 1990 


6. Performing Organization Code 

7. Author(s) 

Robert G. Deissler 

8. Performing Organization Report No. 

E-5568 

10. Work Unit No. 
505-90-01 

9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 

14. Sponsoring Agency Code 

1 

15. Supplementary Notes 

16. Abstract 

The variation of the velocity-derivative skewness of a Navier-Stokes flow as the Reynolds number goes toward 
zero is calculated numerically. The value of the skewness, which has been somewhat controversial, is shown to 
become small at low Reynolds numbers. 

17. Key Words (Suggested by Author(s)) 

Skewness factor; Numerical solution; 
Low-Rev nolds-number; Decaying flow; 
Navier-Stokes; Turbulence 

18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 34 

— i — — 1 

19. Security Classif. (of this report) 20. Security Classif. (of this page) 

Unclassified Unclassified 

21. No. of pages 22. Price' 

10 A02 


NASA FORM i«m oct 86 .p or sa | e by the National Technical Information Service, Springfield, Virginia 22161 






National Aeronautics and fourth CLASS MAN. 

Space Administration 

Lewis Research Center address correction requested 

Cleveland. Ohio 44135 

OMcMBu*wm 

PanaPy br Private Uaa S300 



Postage and Faes Paid 
National Aeronautics and 
Space Administration 
NASA-451 





